Fragmentation Phase Transition in Atomic Clusters IV 
— Liquid— gas transition in finite metal clusters 
and in the bulk — 

D.H.E. Gross and M.E.Madjet 
Hahn-Meitner- Institute Berlin, Bereich Theoretische Physik, Glienickerstr. 100 

14109 Berlin, Germany 
and 

Freie Universitat Berlin, Fachbereich Physik 
Q\ ■ February 1, 2008 



in 

(N 



Within the micro-canonical ensemble it is well possible to identify phase- 
transitions in small systems. The consequences for the understanding of phase 
transitions in general are discussed by studying three realistic examples. 

We present micro-canonical calculations of the fragmentation phase tran- 
sition in Na-, K-, and Fe- clusters of N = 200 to 3000 atoms at a constant 

> ■ 

£^ ■ pressure of 1 atm. The transition is clearly of first order with a back-bending 

micro-canonical caloric curve Tp(E, V(E, P)) = dS(E, V(E, P))/dE\p. From 
j^. \ the Maxwell construction of f3p{E /N ,P) = 1/Tp one can simultaneously 

determine the transition temperature T tr , the specific latent heat qi at , and 
the specific entropy-loss As sur f linked to the creation of intra-phase sur- 
faces. T tr As sur f * NftAirr^gNjff) = 7 gives the surface tension 7. Here 

47rr^ s iVgj^ = J2^i* 47rr^, s m^ //3 is the combined surface area of all fragments 
with a mass rrij > 2 and multiplicity N{. All these characteristic parameters 
are for ~1000 atoms similar to their experimentally known bulk values. This 
finding shows clearly that within micro-canonical thermodynamics phase tran- 
sitions can unambiguously be determined without invoking the thermodynamic 
limit. However, one has carefully to distinguish observables which are defined 
for each phase-space point, like the values of the conserved quantities, from 
thermodynamic quantities like temperature, pressure, chemical potential, and 
also the concept of pure phases, which refer to the volume of the energy shell 
of the N-body phase-space and thus do not refer to a single phase-space point. 

At the same time we present here the first successful microscopic calcula- 
tion of the surface tension in liquid sodium, potassium, and iron at a constant 
pressure of 1 atm.. 

PACS numbers: 05.20.Gg, 05.70.Fh, 64.10.+h, 64.70.Fx, 61.46,+w 



I. INTRODUCTION 



Thermodynamics is commonly believed to address to large many-body systems close to 
the thermodynamic limit. It is further common believe that the three popular ensembles, 
e.g. for simplicity here at given pressure P, the micro-canonical {E, N, P}, the canonical 
{T, N, P}, and the grand-canonical {T,fi,P} agree in that limit and describe the same 
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physical behaviour. But this is not so at first order phase transitions as e.g. the liquid- 
gas transition. Fixing the temperature to the boiling temperature but leaving the energy 
undetermined the canonical ensemble is a mixture of pure liquid and pure gas configurations 
with equal probabilities. Drawing different members of the canonical ensemble, the specific 
heat fluctuates by the amount of the specific latent heat. In the canonical representation one 
would not be able to see a pot of boiling water with a surface dividing the liquid from the 
vapor. This is only possible when the system is confined to a given specific energy e = E/N 
i.e. in the micro-canonical ensemble. The energy controls what portion of the water is liquid 
and what is steam. Every-day macroscopic experience is represented by the micro-canonical 
ensemble because the energy supply is seldomly unrestricted. 

Micro-canonical thermodynamics explores the topology of the N-body phase space and 
determines how its volume 

r E rr fiN-Q 

W{E , N, V) —e s — iEt J/ (Mr _ 6 W*.,*) ~E') (1) 

depends on the fundamental globally conserved quantities namely the total energy E = N*e, 
the angular momentum L, the mass (number of atoms N), charge Z, the linear momentum 
p, and last not least the available spatial volume V of the system. This definition is the 
basic starting point of any thermodynamics since Boltzmann [|l]. It is an entirely mechanistic 
definition. If we do not know more about a complicated interacting N-body system but the 
values of its globally conserved macroscopic parameters the probability to find it in a special 
phase space point (N-body quantum state) is uniform in the accessible phase space. It is 
of course a completely separated and difficult question, outside of thermodynamics, if and 
how a complicated interacting many-body system may explore its entire accessible phase 
space. This question does not concern us here. The present work is an attempt to develop 
the thermodynamics of realistic systems entirely from their mechanics without invoking any 
additional assumption like the use of the thermodynamic limit. 

Before we proceed, we have to emphasise the concept of the statistical ensemble. Each 
phase space cell of size (2irh) 3N ~ 6 /5E corresponds to an individual configuration (event) of 
our system. (The factor 5E is the arbitrarily chosen small width of the energy shell. In 
the following all energy variables, also the temperature, are in units of SE, here in leV.) 
Clearly, the volume e s of the phase space is the sum (ensemble) of all possible phase space 
cells compatible with the values of energy etc.. While the conserved, excessive quantities, 
energy, momentum, number of particles, and charge can be determined for each individual 
configuration of the system, i.e. at each phase-space point or each event, this is not possible 
for the phase space volume e s , i.e. the entropy S(E, V, N) and all its increments like the 
temperature T = (dS(E, V, N)/dE)~ l Q), the pressure P(E, V, N) = TdS(E, V, N)/dV, and 
the chemical potential \i = —TdS(E, V, N)/dN. Their determination demands a measure of 
the phase-space volume e s or its resp. variation. They are ensemble averages. Only in the 
thermodynamic limit, for systems with infinitely many particles N may e.g. the temperature 
be determined in a single configuration by letting the energy flow into a small thermometer. 



1 )This is the thermodynamical definition of temperature. It is in general different from its me- 
chanical one as 2/3 of the mean kinetic energy per particle. This is also true for the pressure. 
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For a finite system, e.g. a finite atomic cluster, the temperature, its entropy, its pressure can 
only be determined as ensemble averages over a large number of individual events. E.g. in 
a fusion of two nuclei one obtains the excitation energy in each event from the ground-state 
Q-values plus the incoming kinetic energy whereas the temperature of the fused compound 
nucleus is determined by measuring the kinetic energy spectrum of decay products which is 
an average over many decays. It is immediately clear that the size of S is a measure of the 
fluctuations of the system. 

In what follows we discuss the most dramatic phenomena in thermodynamics: the occur- 
rence of phase transitions. We will try to interpret them micro-canonically as peculiarities 
of the N-body phase space. We will avoid the concept of the thermodynamic limit as we 
believe that this is not really essential for the understanding of phase transitions. We will see 
that details about the transitions become more transparent in finite systems. Then however, 
one needs a modified definition for the concept of phase transitions. 

In the first paper of this series on micro-canonical thermodynamics and fragmentation 
of atomic clusters (papers I — III, [@-[|]) we introduced a new criterion of phase transitions, 
which avoids any reference to the thermodynamic limit and can also be used for finite sys- 
tems: The anomaly of the micro-canonical caloric equation of state T(E/N) where 
dT/dE < i.e. where the familiar monotonic rise of the temperature with energy is inter- 
rupted. Very early the anomaly of the caloric curve T(E/N) was interpreted as signal for 
a phase transition in small systems in the statistical theory of multi-fragmentation of hot 
nuclei by Gross and collaborators see refs. |5]|| and the review article . Challa and Hether- 
ington came to the same conclusion at about the same time in their paper on the Gaussian 
ensemble [BJ. Within the Gaussian ensemble one can transform smoothly from the micro- 
canonical to the canonical ensemble. They concluded that there is a hierarchy of ensembles 
: micro-canonical — > canonical — > grand-canonical in which the content of information 
about the system diminishes from left to right. They claim it is impossible to recover the 
fundamental micro-canonical ensemble from the canonical. I.e. if the canonical partition 
sum is only known with limited accuracy, may be due to some approximation, then the 
back-transform is dangerously inaccurate. This is a fortiori the case at the most interesting 
points: At phase transitions of first order, see below. Of course, this is because of the high 
sensitivity of the Laplace-back transform from the canonical to the micro-canonical parti- 
tion sum. (E.g. as discussed below, the Laplace transform E — >T suppresses configurations 
with liquid-gas phase separation by the exponentially small factor exp(— N 2 ^a sur f /T tr ). 
The Laplace-back transform must consequently enhance these configurations again by the 
inverse of this factor.) For the case of the nuclear level-density formula by Bethe this was 
demonstrated by ref. ||. Bixton and Jortner [[II]] linked the back-bending of the micro- 
canonical caloric curve to strong bunching in the quantum level structure of the many-body 
system i.e. a sudden opening of new phase space when the energy rises. Their paper offers 
an interesting analytical investigation of this feature. 



2 )It is convenient to consider specific quantities as the energy e = E/N or the specific entropy 
s = S/N per atom. In many cases the micro-canonical caloric curve T(e) depends only weakly on 
the particle number N and the extrapolation to the thermodynamic limit is then easy. 
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FIG. 1. Micro-canonical caloric curve Tp(E/N = e,V(E, P,N),N) at constant pressure (full 
square points), number of fragments Nf r with m; > 2 atoms, and the effective number of surface 
atoms N e jj = J2 m i Ni = total surface area divided by 47rr^, s . The two shaded areas correspond 
to the two equal areas between (3(e) = 1/T(e) and the Maxwell-line (3 = 1/Tt r . In the lower panel 
s(e) = Jq (3(e')de' is shown. In order to make the intruder between E\ and £3 visible, we subtracted 
the linear function 25 + 11. 5e. The double-tangent to s(e) is the concave hull to s(e) between E\ 
and £3. Its derivative is the Maxwell-line in the upper two panels. 

A phase transition of first order is characterised by a sine-like oscillation, a "back- 
bending" of T(e = E/N) c.f. fig.|I]. As shown below, the Maxwell-line which divides the 
oscillation of dS/dE = (3(e) = 1/T into two opposite areas of equal size gives the inverse of 
the transition temperature T tr , its length the specific latent heat qi at , and the shaded area 
under each of the oscillations is the loss of specific entropy As sur f. The latter is connected 
to the creation of macroscopic intra-phase surfaces, which divide mixed configurations into 
separated pieces of different phases, e.g. liquid droplets in the gas or gas bubbles in the 
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liquid. Even nested situations are found in some other cases like liquid droplets inside of 
crystallised pieces which themselves are swimming in the liquid in the case of the solid - 
liquid transition, see e.g. the experiments reported in |jTT|. I.e. at phase transitions of first 
order inhomogeneous "macroscopic or collective" density fluctuations are common, boiling 
water is certainly the best known example. Phase-dividing surfaces of macroscopic size exist 
where many atoms collectively constitute a boundary between two phases which cause the 
reduction of entropy by As sur f. As the entropy is the integral of /3(e): 

s (e) = I* /3(e')de' (2) 
Jo 

it is a concave function of e (d 2 s/de 2 = d/3/de < 0) as long as T(e) = (3~ x shows the usual 
monotonic rise with energy. In the pathological back-bending region of (3(e) the entropy 
s(e) has a convex intruder of depth As sur f P2|(c.f. fig.EIJ;). This has the consequence that 
the specific heat (dT/dE)' 1 becomes negative. This astonishing result was discussed early 
by Thirring ||13|| . At the beginning (> si) of the intruder the specific entropy s(e) is reduced 
compared to its concave hull, which is the tangent to s(e) in the points e± and £3. The 
derivative of the hull to (3(e) follows the Maxwell-line in the interval E\ < e < £3. In the 
middle, (£2), when the separation of the phases is fully established this reduction is maximal 
= As sur f and at the end of the transition (e 3 ) when the intra-phase surface disappears As sur f 
is gained back. Consequently, the two equal areas in (3(e) are the initial loss of surface entropy 



As sur f and the later regain of it. Due to van Hove's theorem |J4| this convex intruder of 
s(e) must disappear for "thermodynamically stable systems" as we consider in this paper 
in the thermodynamic limit which it will do if As sur f ~ iV -1 / 3 . This is why a transition 
of first order may easier be identified in finite systems where the intruder can still be seen. 
The intra-phase surface tension is related to As sur f by 7 sur / = As sur f * N * T tr /surf. -area. 

In conventional thermodynamics a phase transition of first order is indicated by a singu- 
larity (jump) in the specific heat cp(T) in the thermodynamic limit (N — > oo| e= jv/y=const)- 
It is easy to see the relation to our present criterion. The canonical partition sum is obtained 
from the micro-canonical volume W(E) of the accessible N-body phase space, in units of 
(2nh) 3N ~ 6 /5E, or the number of quantum mechanical N-body states, W(E) = e Ns ( £=E / N ) 
by the Laplace transform from the extensive energy to the intensive temperature T = 1/(3: 

Z({3) = / W(E)e~ l3E dE. (3) 



n 



In the thermodynamic limit (N — ► oo| e= jv/i/=const) the partition function Z((3) and the bulk 
specific heat c(T) are obtained via the Laplace transform (D) in saddle-point approximation 
(see e.g. |15| App.2B) : 

POO 

Z((3) = / W(Ne) e~ pm Nde (4) 
Jo 

~ T e N ^-^ ^/27iN/c(e) (5) 
d s(e)\ £ - = {3. (6) 
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The exponent in eq.(|5]) is the free energy —Nf(T)/T. Eq.|6] establishes a relation between 
T and e, the caloric equation of state. Whereas T{e) is single valued everywhere, this is 
not true for its inverse e(T) in the back-bending region. There are in general 3 solutions for 
T = T tr of eq.(|) here c.f. fig.|| e = £\, e 2 , and e 3 . At the second solution e 2 of equation 
(§) s(e) has a positive curvature and the saddle point is a minimum in the direction of 
integration not a maximum. Consequently, the Laplace transform (f|) towards the canonical 
partition function Z(T) jumps over this region which becomes exponentially suppressed 
by a factor e~ A ^ sur f' TtT where Af sur j oc N 2 ^ 3 is the positive surface free energy. In the 
thermodynamic limit (N — > oo) the integration follows the concave hull of s(e) here. Then 
c(T) gets the singularity = qi at 5(T — T tr ) but otherwise remains finite P^B|JT7[] , which is just 
the conventional signal of a phase transition of first order, c.f. Ruelle [y|]. This jump over 
the energies {e\ < e < £3) of configurations with phase-separation with the consequence of 
their exponential suppression is the mathematical reason for the loss of information in the 
canonical ensemble. It is just the crucial information about the phase transition. 

Our approach to phase-transitions of first order is complementary to the conventional 
approaches where the separation of the system into two homogeneous phases by a — in 



general — geometrical interface is investigated, e.g. |T9] . The problems due to the large 



fluctuations of this interface are numerous and severe, see e.g. the discussion on fluid inter- 
faces by Evans [20]. These fluctuation are of course crucial for the interfacial entropy and 
consequently for the surface tension also. The two main differences of our approach are that 

• we do not start with the geometry (planar or spherical) of the interface but focus our 
attention to the entropy of the phase separation. This turns out to be much simpler 
than the geometric approach. Moreover, it is the entropy of the interface that decides 
the transition not its geometrical interpretation. 

• The micro-canonical ensemble allows for large scale spatial inhomogeneities, whereas 
the canonical ensemble suppresses spatial large scale inhomogeneities like phase- 
separations exponentially {in the case of a phase transition of first order oc 
exp(—o~N 2 / 3 /T tr ) , where a is the surface tension parameter (a = 47rr^ s 7) in the liq- 
uid drop parametrisation of the ground-state binding-energy of the clusters here of 
course at the boiling point, r ws is the radius of the Wigner-Seitz cell pT|}. This is an 
example of the information loss in the canonical ensemble mentioned by Challa and 
Hetherington ||. 

This characterisation of phase transitions is purely thermodynamically. We have not 
yet defined what a phase is. Much effort is spent in Ruelle's book to define pure phases 
as those configurations for which in the thermodynamic limit observables survive increasing 
coarse-graining, for which space-averaged quantities do not fluctuate, c.f. chapter 6.5 in 



Ruelle's book ||18|| . Of course, this definition works in the thermodynamic limit only. It does 
not address to finite systems. For a finite system it is not possible to decide if a specific 
configuration corresponds to a pure phase or not. The situation is analogous to the definition 
of the temperature, see above. Again to be a pure phase is a feature of the whole ensemble 
not of a single phase-space point (configuration). We use a statistical definition of a pure 
phase: A configuration belongs to the ensemble of pure phases — including its fluctuations — 
at concave points of S(E, N, V(E, P)) with d 2 S/dE 2 < 0. At the two crossing points of the 
micro-canonical caloric curve Tp(E/N) = T tr , ei, e% we say the system is in the pure "liquid" 
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phase or in the pure "gas" phase respectively. This is at this moment nothing more than a 
working hypothesis. However, we will see that for the systems we have investigated (2-dim 
Potts models, fragmentation of sodium, potassium, and iron clusters with N = 200 — 3000 
atoms) these definitions give very realistic values for all parameters of the transitions quite 
similar to the ones in the bulk. Of course, the parameters depend weakly on the number 
of particles considered. However, there is no reason to emphasise any qualitative difference 
between the bulk transition and the "phase transition" with or without quotation marks in 
finite systems as defined above. 

It is evident, at a phase transition of first order (T = T tr ) the members of a canonical 
ensemble split into two distinct groups with different energy /particle (e\ and £3) but with 
equal probability e~^ T , a "liquid" and a "vapour" phase because both have the same free 
energy. At T tr the energy per particle e fluctuates by the specific latent heat whatever the 
coarse-graining might be. This is illustrated in figfj. Consequently, at a first order phase 
transition the micro- canonical ensemble where energy does not fluctuate is different from the 
canonical one even in the thermodynamic limit at T = T tr .. 

If the shaded area As sur f under the oscillation of (3(e) (see fig.[j]) disappears and if the 
latent heat qi a t = 0, we have a continuous phase transition. In the thermodynamic limit the 
Laplace transform eq. then has only one stationary point at T tr in this case, which is a 
saddle point of the caloric curve T{e) and the specific heat c(T) = de/dT is continuous and 
has a pole at T = T tr , c.f. fig.l of paper I 0. 




FIG. 2. Probability Pt r (s) to find the system in the canonical ensemble at a phase transition of 
first order (T = T tr ) at the energy/particle e. At s = E\ it is " liquid" and at e = £3 it is in the 
pure "gas" phase 

In ref. [|12| we showed for a two-dimensional Potts-model that all the above mentioned 
micro-canonical parameters of a phase transition of first order are within a few percent close 
to their bulk values for relatively small systems like a couple of hundreds spins. This is not 
only true in model systems but also in realistic continuous systems like liquid metals. In 
this paper (IV. of our series) we will follow the micro-canonical fragmentation transition 
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in clusters of sodium, potassium, and iron with increasing number of atoms towards the 
well known liquid-gas transition of the bulk. In contrast to the fragmentation transition 
in isolated clusters which have to be treated at an approximately constant volume defined 
by the short range of dissipation the liquid-gas transition must be studied at constant 
pressure, here at 1 atmosphere. The rapid convergence of the three characteristic parameters, 
Ttr, qiati and o surf towards their bulk values (c.f. tables | and |J) is a further check of the 
validity of our characterisation of a phase transition by the anomaly of the micro-canonical 
caloric curve and of our numerical method (MMMC) to treat the fragmentation of small 
systems like nuclei or atomic clusters. 

The advantage of our new signal of a phase transition is the following: As shown in 



|lq , [T7|Jl^ , |2T|| the micro-canonical caloric equation of state T(E/N) gives this signal of a 
phase transition already for relatively small systems. Moreover, T(E/N) becomes rather 
quickly independent of the number of particles outside the back-bending region, which itself 
degenerates for large N towards the Maxwell-line T = T tr . 
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II. THE MICRO-CANONICAL LIQUID GAS TRANSITION 
AT CONSTANT PRESSURE 



The entropy of the system follows from the fundamental formula 



S{E, V,N) = k* lnW{E, V, N), 



(8) 



where W(E, V, N) is the number of all accessible quantum states of the N-body system in 
the energy interval E,E + 5E and volume V. In what follows we take Boltzmann's constant 
k — 1. The thermodynamic temperature at constant volume is defined by (in the following 
we suppress the dependence of the number of atoms N, as this is hold fixed in all derivatives) 



T(E,V) 



P(E,V) 



dS(E, V) 



dE 



ftr-2 

< — > 

2Ek >J 



(9) 
(10) 



where f tr denotes the number of translational-rotational degrees of freedom and E' is the 
remaining kinetic energy (c.f. eq.B43 pi]). The pressure is given by 



P(E,V)=T(E,V) 



T(E,V) < 



( dS(E,V) ' 

dWmWr 

> 



dV 



w m w r is the spatial weight given in eqs.B25, B30 of [21]: 



N ' 

1 * m .• 



( 1 



6 \27Th 
Nf + Nn 

V 



-i N m N f 



Vrr 



n 



1 4tt 
\ [(2tt^) 3 T 



(R 



s;y.s 



Rjf\/NCC 



> +terms due to the change of NCC, 



(11) 

(12) 
(13) 

(14) 
(15) 



where NCC — 1 is the number of unsuccessful attempts to put the fragments j into the 
given volume. It takes care of the varying avoided volume. In the present calculations at 
a pressure of latm we have NCC ~ 1 so that its variation does not concern us here. At 
higher pressure, more close to the critical point, the pressure due to the change of the avoided 
volume is more important. Nf + N m is the average number of fragments including monomers 
and V is the available volume of the system. (As already discussed above, in contrast to 
here, for the fragmentation phase transition of free finite clusters the calculations were done 
in papers I III P,|i~2|,|4]| at constant freeze-out volume not at constant pressure.) We can 
calculate the thermodynamic temperature at constant pressure by the following: 



T P (E,P) 



Pp(E,V(E,P)) 



dS(E,V(E,P)) d(SP) 



dE 



d(SP) 

d(EV) 

d(EP) ~ °iMn 

V ' d(EV) 



(16) 



Using the Jacobian 

d{SP) 
d(EV) ~ 



dS(E,V) dP(E,V) dS(E,V) dP(E,V) 



dE 



dV 



dV 



dE 



v 



(17) 
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we get 



~dS 


dP 




dS 


dP 




dE 


vdV 


E 


w 


EdE 


V 



dP 




dV 


E 



(18) 



and 



Pp(E, P) = (3(E, V) * 

On the other hand we have from eqs.[Tl]|] : 

dP(E, V) 



P 



dP 




dE 


V 


dP 




dV 


E 



(19) 



BE 
dP(E, V) 



dT(E, V) dS(E, V) 
v 



dE 



dV 



dV 
dT(E, V) 



dT(E, V) dS(E, V) 

E 



dE 



-T 



dV 

2 d(3(E,V) 



dV 



+ T 



+ T 



d 2 S(E, V) 

dEdV 
d 2 S(E,V) 



dE 



v 



-T 



E dV 2 

2 d 2 S(E,V) 



d 2 E 



v 



(20) 
(21) 
(22) 



Using w m and w r given in the appendix B of the review paper [21] and assuming for the 
moment the avoided volume to keep close to the initial cluster volume (at low densities of 
the system), we finally get 



{3 P (E,P) = P(E,V(E,P)) 



T <N t > (< 



Nt(ftr-2) 

IE' 



> -T < N t >< 



(ft. 



-2)(/tr~4) ^ 



(< N t (N t - 1) > -T < {J ^± >< N t >) 
+terms oc dNCC/dV at low e, high P, 



(23) 



where N t is the total number of fragments including neutral monomers. E' is the remaining 
energy and f tr is the number of translational-rotational degrees of freedom. 



One notes that /3(E, V) =< 



(ftr -2) 

2E' 



>. In the case of vaporisation (only monomers), we 



find 



Pp(E,P) = P{E,V(E,P)) 



1 + 



2N 



3N-8 



(24) 



where N is the number of atoms. At very large number of atoms we get the same formula 
as derived for an ideal gas consisting of iV particles 



(3 P (E,P)=(3(E,V(E,P))-. 



(25) 



The corrections oc dNCC/dV take care of the fact that the avoided volume is larger than 
the total eigenvolume of the fragments. It is bigger at low excitation energies when the 
fragments are larger than at higher excitation. These variations, consequently, contribute 
to the pressure (eq.|15|). These corrections turn out to be negligible at 1 atm.. They are, 
however, more important at higher pressure when one approaches the critical point. 

The micro-canonical ensemble with given pressure W(E, P, N) must be distinguished 
from the (in spirit) similar constant pressure ensemble {H,P,N} introduced by Andersen 
|22]-p4 where a molecular-dynamic calculation with the hypothetical Hamiltonian [^ 
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T/2/3 N N N M 

H({ ri , Pi }, V,V) = —-^ i mr i *r i + J2H H^V 1 / 3 ) + —V 2 + P E V (26) 

1=1 i=l j>l 

is suggested. Here V is the volume of the system, taken as an additional explicit degree of 
freedom, {rj,pj} are the coordinates and momenta of the atoms scaled with the factor V 1 / 3 , 
is the intra-atomic two-body potential, and M is a hypothetical mass for the volume 
degree of freedom. P E is the given pressure. The total "energy" H, atoms plus V^-degree of 
freedom, is conserved, not the total energy E of the atoms alone. 

This is different to our micro-canonical approach with given E,V(E, P), N where the 
energy E of the atoms is conserved and the pressure is the correct thermodynamic pressure 
(P(E,V) = T(E,V)dS/dV\E)- At each energy the volume V(E, P) is chosen for all mem- 
bers of the ensemble simultaneously by the condition that T(E,V)dS(E,V)/dV\E of the 
whole ensemble is the correct pressure. In this case there is a unique correlation between 
the energy E and the volume V which does not fluctuate within the ensemble even though 
the pressure is specified. At the given energy this is still the {E,V(E,P,N),N} ensemble. 
Moreover, the entropy cannot be calculated directly by molecular dynamics, whereas this is 
possible in our micro-canonical Metropolis Monte Carlo method. 



III. THE LIQUID-GAS TRANSITION OF SODIUM, POTASSIUM, AND IRON 

The microscopic simulation of the liquid-gas transition in metals is especially difficult. 
Due to the delocalisation of the conductance electrons metals are not bound alone by two- 
body forces but experience long-range many-body interactions. Moreover, at the liquid-gas 
transition the binding changes from metallic to covalent binding. This fact is a main obstacle 
for the conventional treatment by molecular dynamics |p5fl . 

In the macro-micro approach we do not follow each atom like in molecular dynamics, 
the basic particles are the fragments. Their ground-state binding energies are taken from 
experiments. The fragments are spherical and have translational, rotational, and intrinsic 
degrees of freedom. The internal degrees of freedom of the fragments are simulated as pieces 
of bulk matter. The internal density of states is calculated from the internal entropy of the 
fragments. It is taken as the specific bulk entropy s(e) at excitation energies e < e max = e\,oii- 
The bulk entropy can be determined from the experimentally known specific heat of the 



solid/liquid bulk matter |2H]. Ebon is the specific energy where the boiling of bulk matter 



starts. This approximation allows to take important enharmonicities of the internal degrees 



of freedom into account e.g. near to the melting point. Details are discussed in . Then the 
metallic binding poses no difficulty for us and the metal — nonmetal transition is controlled 
in our approach by the increasing fragmentation of the system. This leads to a decreasing 
mean coordination number when the transition is approached from the liquid side while 
the distance to the nearest neighbour remains about the same. Exactly this behaviour was 



recently observed experimentally p7|,p8| |. By using the micro-canonical ensemble we do not 



pre-specify the intra-phase surface and allow it to take any form. Also any fragmentation 
of the interface is allowed. It is the entropy alone which determines the fluctuations of the 
interface. Here we present the first successful microscopic calculation of the surface tension 
in liquid sodium, potassium, and iron. 
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The figure [1] shows the micro-canonical caloric curve T(e) for a system of N = 200 and 
1000 sodium atoms. The back-bending of Tp(e) can be clearly seen. At e ~ e<i ~ 0.7 — 0.8eV 
the thermodynamic temperature drops suddenly due to the rapid increase of the number 
of fragments. This induces a jump in the increment ds/de = 1/T(e) of the entropy s(e). 
The number of fragments with mass rrii > 2 increases slowly up to 62 and from there on 
jumps up and decays continuously down to 0. Arcr^^^ = 4vrr^ s m^ 3 * Ni is the total 
surface area of the fragments. For e < £2 the size rrii of the fragments decreases due to 
an increasing evaporation of monomers, c.f. fig.[|, but the number of fragments increases 
such that the total surface area decreases more weakly with rising excitation. For N = 200 
and 1000 it keeps even approximately constant « 47rr^ s iV 2//3 . A more detailed investigation 
shows that in sodium as well as also in potassium we actually may have two interfering 
transitions: One from evaporation of monomers and smaller fragments with a large residue 
towards multi-fragmentation with several medium sized fragments at ~ 0.7eV/atom and a 
second one from multi-fragmentation into a pure gas of monomers at e > 1.2eV/atom. For 
systems with more atoms, iVo > 3000 c.f. fig.|3], the multi-fragmentation in sodium moves 
towards larger excitation and melts together with the vaporisation dip at ~ l.l-1.2eV/atom. 
In the canonical ensemble all these important details become suppressed or even hidden. 
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FIG. 3. Same as |T| but for Na^Q 00 . The four small figures at the top show the mass distribution 
of fragments at four different excitation energies which are indicated in the main figure by their 
number. The small vertical numbers on top of the mass-distributions give the real number of 
fragments e.g.: 2:7.876 means there are 7.876 dimers on average at e = 0.442eV/atom. 



The bulk values of a are calculated from the experimental surface tension 7 by: 



O'bulk 



ws 1 

w' 3vM 



30 ' 



4vrL 



a bulk (T) = 6.242 * 10~ 5 * ^{r 2 wsl )\ TmM 1 



1 d-) 
^dT 



2 dv 
3v dT 



(T — T, 



melt ) 



(27) 
(28) 

(29) 



where r ws in [A] is the Wigner-Seitz radius at melting, v in [^-] is the specific volume at 
melting, L is the number atoms per kg- mole (Loschmid's number), and M is the molecular 
weight. 7 | T 



dj / 'dT\T meH are the experimental surface tension in [^P-] and v |t w< ) Jf|t,, 



T-elf dTl 1 melt 



are the specific volume in [|— ] and its temperature derivative in [ at the melting point 

given by |29| . The values by Iida and Gutherie [29] are slightly different from the values from 
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Miedema | 30| . The uncertainty of the experimental values for the bulk surface tension at 
the melting point are quoted |29 to be ~ 5 — 10% at T m whereas the extrapolation towards 
T = resp. towards the boiling point is done with the parameters -^\ me it and ~^\ me it 
which have an estimated error of ±50%. We give the in this way extrapolated values of 
(Jbuik{Tbott) for iron in the last column of table [TI[ 

Inspecting these numbers we find the liquid-drop parameter a s = er(T = 0) = 47rr£, s * 
7|t=o to be less by about 30% than the values of a s determined for the ground-state binding 
energies of real clusters averaged over the shell effects in |}1 which we have of course used 
in our calculation for the ground-state binding energies of the fragments. The origin of this 
experimental discrepancies is yet unknown [[31]]. Consequently, we think we should compare 
our theoretical values for ab a with values for GbuikiT) which are consistent at T = with 
the values given for a s by |31] and than extrapolated to T = T^oii with formula ^9|. These 
are listed in the last column of table PI 
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Excitation Energy (eV/ atom) 

FIG. 4. Micro-canonical caloric curve Tp(E/N = e, P) (full square points), number of monomers 
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To check this conclusion we recalculated the transition parameters for Fei~ 000 and Fe^Qo 
using instead of a s = 5.1eV, a s = 4.0eV , the value of (Jbuik(T = 0) as estimated with the 



help of formula [29| from the experimental surface tension of the bulk at melting. Again 
the resulting transition parameters listed in TABLE |J approach nicely the corresponding 
"experimental" bulk values at T = Tb u- As our theory uses the ground-state (T = 0) 
binding energies for the clusters as input values, the theory predicts not the total surface 
tension but its temperature dependence, the so called "entropic part of the surface tension" . 
The good reproduction of the adjusted value at boiling in TABLE | as well as the true one 
in TABLE [TJ shows the high quality of our method. 

Fig. [| shows the same caloric curve Tp(e) as fig.[l] but now the number of evaporated 
monomers. At e > 1.2eV the system is totally vaporised into monomers. At £2 ~ 0.7eV the 
character of the decay of the system changes and this can also be seen in N mon (e). 
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FIG. 5. Micro-canonical caloric curve Tp(E/N = e, P) (full square points), number of monomers 



The decay of potassium is in all details similar to that of sodium. Therefore we don't 
show here the corresponding figures. The liquid-gas transition in iron is different from that 
of the alkali metals: Due to the considerably larger surface energy parameter a s in the liquid 
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drop formula of the ground-state binding energies of iron compared to alkali metals there 
is no multi-fragmentation of iron clusters at P = 1 atm. Iron cluster of N < 3000 atoms 
decay by multiple monomer evaporation c.f. fig|| The value of given in tables | and 

|TT] for iron is taken as the average of over the energy interval £1, £2- 

Table | gives a summary of all theoretical parameters for the liquid-gas transition in 
clusters of iVo = 200 — 3000 Na, K, and Fe atoms and compared with their experimental 
bulk values. The transition-temperature T tr , the specific latent heat and the entropy 
gain of an evaporated atom s^oii are well approaching the experimental bulk values. As sur f 
is the area under the back-bend of (3{e). As sur f * N is the total entropy loss due to the 
interfaces equal to 

£ ^r 2 ws m* /3 N a /T tr = N%)a/T tr . (30) 





N 


200 


1000 


3000 


bulk 




T tr [K] 


816 


866 


948 


1156 




Qlat [eV] 


0.791 


0.871 


0.91 


0.923 


Na 


Sboil 


11.25 


11.67 


11.2 


9.267 




As sur f 


0.55 


0.56 


0.45 






iV e/.f 


39.94 


98.53 


186.6 


00 




^/Ttr 


2.75 


5.68 


7.07 


7.41 




Ttr [K] 


697 


767 


832 


1033 




Qlat [eV] 


0.62 


0.7 


0.73 


0.80 


K 


Sboil 


10.35 


10.59 


10.15 


8.99 




As sur f 


0.65 


0.65 


0.38 






iV e./7 


32.52 


92.01 


187 


00 




(y/Ttr 


3.99 


7.06 


6.06 


7.31 




Ttr [K] 


2600 


2910 


2971 


3158 




Qlat [eV] 


2.77 


3.18 


3.34 


3.55 


Fe 


Sboil 


12.38 


12.68 


13.1 


13.04 




As sur f 


0.75 


0.58 


0.77 








22.29 


65.40 


142.12 


00 




^/Ttr 


6.73 


8.87 


16.25 


17.49 



TABLE I. Parameters of the liquid-gas transition at constant pressure of latm. in a mi- 
cro-canonical system of iVo interacting atoms and in the bulk. Sb n = qi a t/Tt r , it is interesting 
that the value of Sboil for all three systems and at all sizes is near to SboU = 10 as proposed by 



the empirical Trouton's rule [Q, N^A = J2 m i /0 ^i an d <y/T tr = NQAs sur f/N^^ 
ues a '/T tT are adjusted according to formula |29] to the input values of a s taken for the T 



T 2/3 



The bulk val- 


liquid-drop surface parameters from ref. [31] which we used in the present calculation for the 
ground-state binding energies of the fragments. 
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N 


1000 


3000 


bulk 




T tr [K] 


2994 


3044 


3158 




qiat [ev\ 


o.lo 


6.1 1 


O.OO 


Fe 


Sboil 


12.13 


12.47 


13.04 






0.48 


0.39 






N 2/ ' S 


65.74 


136.78 


oo 




cr/Ttr 


7.30 


8.62 


8.86 



TABLE II. Same as table |, but using for the ground-state liquid-drop parameter a s the value 
compiled by extrapolating the experimental surface tension of bulk iron at melting down to T = 
(c.f. the discussion in the text) or to T = T\, n. 
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Of course, the transition temperature T tr and the latent heat qi at of small clusters are smaller 
than the bulk values because the average coordination number of an atom at the surface of 
a small cluster is smaller than at a planar surface of the bulk. It is somewhat surprising that 
the surface tension is rising from N Q = 200 to N = 3000. One might think it should fall as 
a is expected to be proportional to the number of broken bonds per surface area. However, 
this is a purely static argument and does not take the entropic (ffuctuational) part of a into 
account. 

IV. THE RELATION TO THE METHOD OF THE GIBBS-ENSEMBLE 

The alternative method to simulate the liquid-gas transition within the canonical en- 
semble is by the method of the Gibbs-ensemble [|19| . As configurations with two coexistent 
phases separated by an interface become asymptotically suppressed, see above, in the Gibbs- 
ensemble method the system is sampled in two separated containers which can communicate 
via particle and energy exchange. In one container the system is close to the bulk liquid 
density and specific energy and in the other it is in the corresponding gas phase without 
establishing an inner, and possibly fluctuating interface. Then the transition temperature 
and latent heat can be determined. 

The Gibbs-ensemble method is used to investigate the critical behaviour of systems 
interacting via short range forces like Lennard- Jones potentials, c.f. the recent article by 



The surface tension cannot be determined by the Gibbs-ensemble method. The dramatic 
fluctuations of the physical interfaces between both phases at T = T tr which can be seen e.g. 
in boiling water are a clear warning of the much more complicated and strongly fluctuating 
microscopic topology of the interfaces. They are considerably more complicated than the 
geometry assumed in the Gibbs-ensemble method. The surface tension is proportional to 
the entropy of the interface-fluctuations and consequently it depends essentially on these 
fluctuations. 

V. ALTERNATIVE MICROSCOPIC METHODS TO CALCULATE THE 

SURFACE TENSION 

Microscopic canonical calculations of the surface tension have been done for Lennard- 
Jones systems. A detailed survey is given by Salomons and Mareschal [0. They discuss 
the following alternatives to determine the surface tension: 

1 . from the virial expansion of the free-energy as function of volume conserving distortions 
of the interface: 



where r the interatomic coordinates of the atoms, A is the surface area and 

4>'{fij) is the derivative of the interatomic two-body potential, 



Bruce @. 




(31) 
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2. from histograms of the free-energy as function of potential energy changes due to 
volume conserving deformations. 

The authors of also determine the "entropic" part of the surface tension dj/dT either 
directly from a compilation of j(T) at various temperatures or via the fluctuation equation 

^ = ^ <{ U- Uh> , (32) 

where U is the total potential energy and U its mean value. 

Evidently, the virial expansion works only for systems with two-body forces. So we 
cannot use it for metallic systems. The second method has not been tried either. 



VI. CRITICISM AND NECESSARY IMPROVEMENTS OF THE 
COMPUTATIONAL METHOD 



Several simplifications were made: 

• We use classical not quantum statistics. - - We do not believe that the problems we 
discussed in the present paper depend essentially on that difference. 

• The fragments were assumed to be spherical. — This simplification means that surface 
degrees of freedom of the fragments are treated like bulk excitations. The size of 
the fragments was assumed to be independent of the excitation energy at e < e max . 
Moreover, the static deformation of dimers and trimers is quite important for decay 
thresholds. An improvement of these approximations was not possible in the present 
calculation. 

• Only internal excitations of the fragments of e < e max were allowed. — This approx- 



imation may be justified for the decay of hot finite clusters |26| , |21| where we argued 
that only a transient equilibrium of the decaying cluster over a time r equ is achieved 
and states which live considerably shorter should not contribute as independent decay 
channels. This argument does not work in the present scenario. However, we believe 
that highly unbound states above the boiling excitation e^i in small fragments are 
anyway better described by assuming an independent motion of the decay products, 
just as we do it here. Doubling the actual value of e max had in most cases only little 
influence on the transition parameters obtained. 

The calculation of the total intra-phase surface area is problematic. Consequently, the 
determination of the surface tension per surface area 7 suffers under a great uncer- 
tainty, whereas the determination of the surface entropy As sur f is easy. The situation 



is complementary to the Gibbs ensemble method ||19|| . However, the surface entropy is 
crucial for the occurrence and classification of a phase transition, not its geometrical 
interpretation. 

The input values of the liquid drop parameters of the ground-state binding energies 
of the fragments given by |TJJ suffer from the uncertainty of the temperature of the 



clusters where the experimental data were taken. The bulk measurements of the 



surface tension suffers from many effects c.f. Pq . f29|1 : 
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— The surface tension of alkali clusters might be reduced by oxygen contamination. 

— The volume expansion coefficient ^dv/dT is only known within ±50%. 

— j(T = 0) refers to the solid phase which is most likely anisotropic and cannot 
easily be extrapolated from the melting point. 



VII. CONCLUSION 



The liquid-gas transition in metals at normal pressure of 1 atm. is experimentally well 
explored. Therefore, it is a good test case for our new ideas and computational methods 
of micro-canonical thermodynamics. By allowing a system of N = 200 — 3000 atoms to 
condense or fragment into an arbitrary number of spherical fragment clusters which are 
internally excited to energies of e < e^i per atom and into an arbitrary number of free 
atoms under a fixed given external pressure of 1 atmosphere we directly calculated the 
micro-canonical caloric curve (3(e) = ds/de =< d/de > (c.f. the detailed formula eq . (|23|) ) 
by micro-canonical Monte Carlo methods iMMMC) |^TJ. The volume W(E) of the total 
accessible N-body phase space can then be determined from (3p(e) by s(e) = J £ (3(e')de' and 
W(E) = exp[Ns(E/N)]. The anomaly of the micro-canonical caloric curve Tp(e) signals 
the liquid-gas phase-transition in the finite system. The characteristic parameters T tr , qi at , 
and the surface tension a sur f approach already for ~ 1000 — 3000 atoms the experimentally 
known values of the bulk liquid-gas transition. 

This result is remarkable for several reasons: 

1. It proves that a phase transition of first order in a realistic continuous system can very 
well be seen and classified in small mesoscopic clusters without invoking the thermo- 
dynamic limit. In fact, with the non-vanishing back-bending of Tp(E/N, V(E, N, P)) 
the transition is easier recognisable than in the thermodynamic limit. For the 10-states 
Potts model we know that the micro-canonical parameters depend less on the size of 



the system than the canonical ones ||12|| . Finite size scaling for the micro-canonical 
transition parameters are not known in general. However, here we want to stress that 
there is no reason at all not to discuss phase transitions in small systems. 

Even though the values determined for the specific surface tension a = 47rr^ s 7 are 
somewhat uncertain because of the difficulty to fix and determine the surface area, the 
surface entropy NAs sur f which is the relevant quantity that determines the nature of 
the transition can very well be calculated. 

Even for a realistic metallic system with its long-range many-body interactions the 
Mcrocanoncal Metropolis Monte Carlo simulation method iMMMC) is able to describe 
the liquid-gas phase transitions quite well. This is possible because we do not use 
molecular dynamics with a two-body Hamiltonian but use the experimental ground- 
state binding energies of the fragment clusters which of course take care of the metallic 
bonding of their constituents. 

The intra-phase surface entropy can be microscopically calculated. (Essentially we 
calculated the temperature dependent part of it. The value at T = is an input pa- 
rameter of our calculation.) The surface tension per surface area can be determined if 
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the intra-phase area is known. For the surface entropy the fluctuation and fragmenta- 
tion of the surfaces are essential. At the liquid-gas transition of sodium and potassium 
clusters of sizes as considered here at normal pressure the intra-phase fluctuations are 
mainly due to strong inhomogeneities and clusterization and not due to homogeneous 
stretching. This is consistent with recent experimental evidence pTp8 |. 



5. Micro-canonical thermodynamics gives an important insight into the details of the 
transition which are hidden in the canonical approach. The transition in sodium is a 
good example: Figs.[3] and |^ show how the transition is a combination of actually two 
transitions, one from evaporation of very light fragments out of a large residue into a 
multi-fragmentation of the system into many larger pieces, and second the complete 
decay into a gas of monomers. 

6. The success of the Mcrocanoncal Metropolis Monte Carlo sampling method to repro- 
duce the known infinite matter values of the liquid-gas transition is also a promising 
and necessary test of MMMC to describe nuclear fragmentation correctly [H,PTf and 
to be able to get insight into the (eventually critical) behaviour of the nuclear matter 
liquid-gas transition. This is important as an experimental test of the predictions of 
the model for nuclear matter is not available. 

7. The division of the macroscopic time-independent observables of the N-body system 
into 

(a) observables which can be determined at each phase-space point, in each individual 
realisation of the micro-canonical ensemble: The globally conserved quantities 
like the energy, the number of particles, the charge, the momentum, the angular 
momentum, and 

(b) into the thermodynamic observables which refer to the size of the ensemble, the 
volume e s of the energy shell of the N-body phase space, and which cannot be 
determined at a single phase-space point, in a single event, like the entropy, the 
temperature, the pressure, the chemical potential 

is very essential. Also the concept of a pure phase and of a phase-transition belongs 
to the second group. 

8. Last not least, there is a very fundamental difference between micro-canonical and 
canonical ensembles. This cannot be emphasised enough: As said above, if one wants 
to study a configuration where two-phases coexist simultaneously and are separated 
by an interface one has to control energy (if pressure is fixed). Thus the the micro 
ensemble {E,P,N] is the adequate ensemble at phase transitions of first order in the 
bulk, not the canonical, at least if the energy resources are not unlimited as they usually 
are. 



A preliminary version of this paper was already published in [BH . 
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